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Abstract 
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1 Polymerized membranes 

Mmembranes have attracted much interest during the last years, especially due to their relevance for biological 
systems. Most of the membranes encountered there are fluid. In this lecture I shall concentrate on another 
fascinating class of membranes, namely polymerized tethered also called solid membranes. These membrane 
have a fixed and constant internal metric. They are realized in experiments (e.g. the spectrin network of red blood 
cells [1,2], or sheets of graphite oxide [3,4]). 

A microscopic model is given by the so-called "spring and bead model" which consists of balls (beads) which 
are connected by springs and form a regular lattice. The difficulty is to incororate the self-avoidance of the beads, 
i.e. the fact that beads cannot intersect. 

2 Representation as non-local field-theory 

We start from the continuous model for a D-dimensional flexible polymerized membrane introduced in [5,6], and 
further studied in [25-28]. This model is a simple extension of the well known Edwards' model for continuous 
chains. The membrane fluctuates in d-dimensional space. Points in the membrane are labeled by coordinates 
X £ and the configuration of the membrane in physical space is described by the field r : x G MP — > 
r(x) £ R"*, i.e. from now on we note r instead of r. At high temperatures the free energy for a configuration is 
given by the (properly scaled) Hamiltonian 

The integral runs over _D-dimensional space and V is the usual gradient operator. The normalizations are 

■- d°x with Sd = and 5'^{r{x) - r{y)) = {4TTf^'^5'^{r{x) - r{y)). The latter term is 

normally used in Fourier-representation S'''(r{x) — r{y)) = e'*'''"'-^-'"'''''", where the normalization of is 

given by = tt"'*''^ / A'^p to have JpS~^ " = cl~'''^'^ . All normalizations are chosen in order to simplify 
the calculations. The first term in the Hamiltonian is a Gaussian elastic energy which is known to describe the 
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Figure 2.1: The critical curve e{D, d) = 0. The dashed line corresponds to the standard polymer perturbation theory, critical in d = 4. 



free "phantom" surface. The interaction term corresponds (for bo > 0) to a weak repulsive interaction upon 
contact. The expectation values of physical observables are obtained by performing the average over all field- 
configurations r(x) with the Boltzmann weight e~^^^\ This average can not be calculated exactly, but one can 
expand about the configurations of a phantom, i.e. non-interacting surface. 

Such a perturbation theory is constructed by performing the series expansion in powers of the coupling con- 
stant b. This expansion suffers from ultraviolet (UV) divergences which have to be removed by renormalization 
and which are treated by dimensional regularization, i.e. analytical continuation in D and d. Long-range infrared 
(IR) divergences also appear. They can be eliminated by using a finite membrane, or by studying translation- 
ally invariant observables, whose perturbative expansion is also IR-finite in the thermodynamic limit (infinite 
membrane). Such observables are "neutral" products of vertex operators 

JV JV 

0=Y[ e*'=<''-(^°) , ^ fca = . (2.2) 

a=l a=l 

Let US now analyze the theory by power-counting. We use internal units /i ~ 1/x, and note [x]^ = 1, and 
[n]^ = — [n]^ = —1 . The dimension of the field and of the coupling-constant are: 

'^:=K = ^^. e:=[bo]^ = 2D-,.d. (2.3) 

The interaction is relevant for e > 0, see figure 2.1. Perturbation theory is then expected to be UV-finite except 
for subtractions associated to relevant operators. We shall come back to this point later. 

For clarity, we represent graphically the different interaction terms which have to be considered. The local 
operators are 

1 = 1, i(vr(x))'=^. (2.4) 

The bi-local operator, the dipole, is 

d''[r{x)-r{y)) = . .. (2.5) 

The expectation-value of an observable is 

{0[r])^ = Norm j V[r] 0[r] e"^M . (2.6) 
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with Norm chosen such that (1);, = 1. Perturbatively, all expectation-values are taken with respect to the free 
theory (again normalized s.t. = 1): 

{0[r])^ = Nona J V[r] 0[r] Q-^hyh 5(Vr-(x))= (2.7) 

A typical term in the expansion of (2.6) is 

{~bZ,^^n" II ...Jl {O. ..... (2.8) 

where the integral runs over the positions of all dipole-endpoints. 

3 Locality of Divergences and the multllocal operator product 

expansion (l\/IOPE) 

The next step to show ist that all divergences are short distance divergences. Note that even for massless theories 
and in the absence of IR-divergences, this is not trivial. Divergences could as well appear, when some of the 
distances involved become equal, or multiple of each other. A simple counter-example is the integral of | |o| — 

\b\\ , where a and b are two of the distances involved. Due to a lack of space, we can not give the proof here. 
The interested reader should consult [7]. 

As in local field-theoies, divergences can be analyzed via an operator product expansion. Intuitively, in the 
context of multilocal theories - by which we mean that the interaction depends on more than one point - we 

also expect multilocal operators to appear in such an operator product expansion, which therefore will be called 
"multi-local operator product expansion" (MOPE) [8,9]. An example for such an operator product expansion is 




(3.1) 

Let us explain the formula. We consider n dipoles (here n = 5) and we separate the 2n end-points into m subsets 
(here m = 3) delimited by the dashed lines. The MOPE describes how the product of these n dipoles behaves 
when the points inside each of the m subsets are contracted towards a single point zj . The result is a sum over 
multilocal operators ^i{zi, . . . , Zm), depending on the rn points zi, . . . , 2™, of the form 

"^dixi - Zl,. . .)^i{zi,Z2,. ■■ ,Zm) , (3.2) 

where the MOPE-coefficients Ci{xi — zi, . . .) depend only on the distances xi — Zj inside each subset. This 
expansion is valid as an operator-identity, i.e. inserted into any expectation value and in the limit of small dis- 
tances between contracted points. No other operator should appear at the points towards which the 
operators are contracted. As the Hamiltonian (2.1) does not contain a mass-scale, the MOPE-coefficients are 
homogeneous functions of the relative positions between the contracted points, with the degree of homogeneity 
given by dimensional analysis. In the case considered here, where n dipoles are contracted to an operator this 
degree is —rwd — [^i\x- This means that 

C^{\{xi - 2i), . . .) = A-"^'*-[*'i-C,(x-i - 21, . . .) , (3.3) 

where ^i\x is the canonical dimension of the operator and — d(2 — /))/2 is the canonical dimension of the 
dipole. 

In order to evaluate the associated singularity, one finally has to integrate over all relative distances inside 
each subset. This gives an additional scale factor with degree D(2n — m). A singular configuration, such as in 
Eq. (3.1), will be UV-divergent if this degree of divergence 

D{2n - m) - n'^-^d - , (3.4) 

is negative. It is superficially divergent if the degree is zero and convergent otherwise. The idea of renormalization 
is to remove exactly these superficially divergent contributions recursively. 
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4 Evaluation of the MOPE-coefficients 



The MOPE gives a convenient and powerful tool to calculate the dominant and all subdominant contributions 
from singular configurations. In this section, we explain how to calculate the MOPE-coefflcients on some explicit 
examples. These examples will turn out to be the necessary diagrams at 1-loop order. 

In the following we shall use the notion of normal-ordering. The first thing, which we use, is that 

.^ikr(x), _ ^ikr(x) ^^^^ 

Explicitly, tadpole-like contributions which are powers of J d^p -p are omitted. This is done via a finite part pre- 
scription (analytic continuation, dimensional regularization), valid for infinite membranes, for which the normal- 
order prescription is defined. Let us stress that this is a pure technical trick, which can be cirumvented at the 
expense of more cumbersome calculations. 

The key-formula for all further manipulations is 

. ^ikr{x) , , ^ipriy) , _ ^kpC(x -y) , ^ikr(x) ^ipr(y) , (4 2) 

This can be proven as follows: Consider the (free) expectation value of any observable O times the operators 
of Eq. (4.2). Then the left- and right-hand sides of the above equation read C = (o le'*''*"^) : :e*f''(«) : and 

n = e'='"='(^-«^ (O :e*'='-(^)e*'"'(^' First of all, for O = 1, the desired equaUty of C = TZ holds, because 
^.^ikr(x)^ipr(y) .^^ ^ ^ ^.^ikr{x) . .^ipr(y) .^^ ^ ^kpC(x-y) _ considcr a non-trivial observable O, and 

contract all its fields r with e'*"'^'"' or e''"'^^', before contracting any of the fields r{x) with r{y). The result is a 
product of correlation-functions between the points in O and x or y, and these are equivalent for both £ and TZ. 
However, contracting an arbitrary number of times e'*'''^^\ leaves the exponential e'*'''^'"^ invariant. Completing 
the contractions for £ therefore yields a factor of c'=p'^(^-y\ and the latter one also appears in TZ. Thus, the 
equality of £. and TZ holds for all O and this proves Eq. (4.2). 

Now proceed to the first explicit example, the contraction of a single dipole with endpoints x and y: 

,^ikr{x) _^-ikr(y) _ 

This configuration may have divergences when x and y come close together. Let us stress that in contrast to <^''- 
theory, these divergences are not obtained as a finite sum of products of correlators: Since C(a;—j/) = |a;— yp"^, 
the latter is always well-behaved at x = y. The singularity only appears when summing an infinite series of 
diagrams as we will do now. To this purpose, we first normal-order the two exponentials using Eq. (4.2) 



.^ik[r{x)-r(y)]_ ^-k \x-y\ " ^) 



Note that the operators e''^'"'^-' and e are free of divergences upon approaching each other, since no more 

contractions can be made. The divergence is captured in the factor e~* l'^"^' . Therefore, we can expand the 
exponential :e''^[''(^'~''(^^l : for small x — y and consequently in powers of \r{x) — r{y)]. This expansion is 

^ : |l + ifc [r{x) - r{y)] - ^ (fe [r{x) - r{y)]f + . . .| : e-^'l^^-^l'" . (4.5) 

We truncated the expansion after the third term. It will turn out later that this is sufficient, since subsequent 
terms in the expansion are proportional to irrelevant operators for which the integral over the MOPE-coefflcient 
is UV-convergent. 

Due to the symmetry of the integration over k the term linear in k vanishes. Also due to symmetry, the next 
term can be simpUfied with the result 



^~ 2d ' ^^^^^ ~ ^^^''^ ■ ■ ■ 



e"* l"'-^! " . (4.6) 



Finally, the integration over k can be performed. Recall that normalizations were chosen such that /j. e 
s'"^^^ to obtain 

"-^''1).+ 



+/3 + . . . , (4.7) 
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where we used the notation = ^{dar){d/3r) and the MOPE-coefficients (reminding Feynman's bra and 

ket notation) 

Ix-vr" (4.8) 

(4.9) 



a^/jj = -^{x -y)c{x -y)p\x -y\ ''''^+^^ 



As long as the angular average is taken (and this will be the case when integrating the MOPE-coefflcient to obtain 
the divergence), we can replace in Eq. (4.7) a-^-is by -f := i(Vr)^ and Eq. (4.9) by 

-i-^=-^\x-yf-'''. (4.10) 

Next consider a real multi-local example of an operator-product expansion, namely the contraction of two 
dipoles towards a single dipole: 



x+u/2 ; 7» , y+v/2 
x—u/2 ^^^^ y — v/2 



I ^ik\r{x+u/2)-r(y+v/2)] f ^ip[r{x-u/2)-r(y-v/2)] (4 11) 

Jk Jp 



This has to be analyzed for small u and v, in order to control the divergences in the latter distances. As above, we 
normal-order operators which are approached, yielding 

^ikr(x-{-u/2) ^ipr(x—u/2) ^^ikr(x+u/2) ^ ^^ipr(x — u/2) ^ 

_ _^ikr(x+u/2)^ipr(x-u/2) _ ^kpC(u) (4 12) 

A similar formula holds when approaching e-''='-(!'+"/2) and e-*^''*^-''/^) . Eq. (4. 1 1) then becomes 

,^ikr(_x+u/2)+ipr(x—u/2)^ _ ^—ikr(y+v /2)—ipr{y—v /2) _ ^kp[C{u)+C{v)] (4 13) 



k J p 

In order to keep things as simple as possible, let us first extract the leading contribution before analyzing sublead- 
ing corrections. This leading contribution is obtained when expanding the exponential operators (here exemplified 
for the second one) as 

.^-ikr(y+v/2)^-ipr(y-v/2) _ _ .^-i(k+p)r(y) (X-|-0(Vr)): (4.14) 

and dropping terms of order Vr. This simpUfies Eq. (4.13) to 

.^i(k+p)r(x), .^-i(k+p)r(y) , ^kplC(u)+C(v)] _ ^^ -^^y 



k J p 



In the next step, first k and second p are shifted k — > k — p then p — > P + f • The result is (dropping the 
normal-ordering according to Eq. (4.1)) 

f ^tk[r(x)-r{y)] f(lk^-p'^)[C(u)+C(v)]_ 
Jk Jp 

The factor of /j. e'*'''"'^'"'''^'! is again a 5-distribution, and the leading term of the short distance expansion of 
Eq. (4.16). Derivatives of the (5-distribution appear when expanding e^2'=^-J'^)[C'(")+'^('')l in k^; these are less 
relevant and only the first sub-leading term will be displayed for illustration: 

j ^^klr(x)-r(y)] f ^-p^ (.u) + C(v)] U+t^ [C (u) + C {v)] + . 

= (^0)|") ~ + tO)|— ) (4.17) 
where in analogy to Eqs. (4.7) to (4.9) 



= [C{u) + C{v)]-'"' , 
= i [^(«) + Civ)]'-"^' (4.18) 
with . — . = 5'^{r{x) - r{y)) and — = {-Ar)S'^{r{x) - r{y)). 
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Let us already mention that the leading contribution proportional to the 5-distribution will renormalize the 
coupling-constant, and that the next-to-leading term is irrelevant and can be neglected. The same holds true for 
the additional term proportional to ( Vr) which was dropped in Eq. (4. 14). 

There is one more possible divergent contribution at the 1-loop level, namely . We now show that the 

leading term of its expansion, which is expected to be proportional to • — • , is trivial. To this aim consider 

P^Ciy-z) ^kplCix-z)-Cix-y)] 

We want to study the contraction of x, y, and z, and look for all contributions which are proportional to • • = 

: e''"'^"^ : : Q-^''<i^+y+^)/3) xhe key-observation is that in Eq. (4.19) the leading term is obtained by 
approximating e'^^I'^'^'^"^)"'^'^"'"^^! w 1. All subsequent terms yield factors of k, which after integration over k 
give derivatives of the ^''-distribution. The result is that 

UOH-(Q 

This means that divergences of . j^-^ are already taken into account by a proper treatment of the divergences 
in Q , analyzed in Eq. (4.7). 



X y z 

=/ 

•/ k.p 



_^ifcr('u), ,Q—ikr{x)^ .^Wfiu). 



k,p 



(4.19) 



5 Renormalization at 1-loop order 

In the last sections, we discussed how divergences occur, how their general structure is obtained by the MOPE, and 
how the MOPE-coefficients are calculated. In the next step, the theory shall be renormaUzed. The basic idea is to 
identify the divergences through the MOPE, and then to introduce counter-terms which subtract these divergences. 
These counter-terms are nothing else than integrals over the MOPE-cocfficients, properly regularized, i.e. cut off. 
We introcuce two renormalization group factors Z (renormalizing the field r) and Zb (renormalizing the coupling 
b). Recalling Eq. (2.1), this is 

= ^ J^K^rix))' +bZ,n^ JJj'[rix)-riy)) , (5.1) 

where r and h are the renormaUzed field and renormalized dimensionless coupling constant, and /i — is the 
renormalization momentum scale. (To be precise: The field r in Eq. (2.1) is the bare field and should be noted 
ro.) 

Let us start to eliminate the divergences in the case, where the end-points {x, y) of a single dipole are con- 
tracted towards a point (taken here to be the center-of-mass z = {x + y)/2). The MOPE was given in Eq. (4.7) ff 
We there have to distinguish between counter-terms for relevant operators and those for marginal operators. The 
former can be defined by analytic continuation, while the latter require a subtraction scale. Indeed, the divergence 
proportional to 1 is given by the integral 



iA-l<|x-j,|<i V^^'^^'f / A-1 x D-s\ ) 



(5.2) 



where A is a high-momentum UV-regulator and L a large distance regulator. For e ~ this is UV-divergent but 
IR-convergent. The simplest way to subtract this divergence is therefore to replace the dipole operator by 

• • • • - • • , (5.3) 



where ^ * • = |a; — j/| This amounts to adding to the bare Hamiltonian (2.1) the UV-divergent counter- 
term 

Afii =-hZkn' j j \x-y\-"\ (5.4) 

J X J y 



which is a pure number and thus does not change the expectation-value of any physical observable. 

We next consider marginal operators: In the MOPE of Eq. (4.9), the integral over the relative distance of 

f ( U+s ) a+B is logarithmically divergent at £ = 0. In order to find the appropriate counter-term, we 
use dimensional regularization, i.e. set e > 0. An IR-cutoff L, or equivalently a subtraction momentum scale 
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H = L ^, has to be introduced in order to define the subtraction operation. As a general rule, let us integrate 
over all distances appearing in the MOPE-coefficient, bounded by the subtraction scale L = Defining 



/ 

J\x 



IL J\x-y\<L 

we need the following counter-term in the Hamiltonian 



a-f /3 



(5.5) 



(5.6) 



subtracting explicitly the divergence in the integrals. The reader is invited to verify this explicitly on the example 
of the expectation value of O = e'''''"^'''"'''^*'!. The solution is given in appendix H of [7]. 

Since the angular integration in Eq. (5.5) reduces to ^ , we can replace Eq. (5.6) by the equivalent 

expression 



for which the numerical value of the diagram is calculated as 



L 



x — y\<iL 



1 

'2D 



Ax 2D-vd 
X 

X 



1 U 
'2D~ 



(5.7) 



(5.8) 



We can now subtract this term in a minimal subtraction scheme (MS). The internal dimension of the membrane 
D is kept fixed and (5.8) is expanded as a Laurent series in e, which here starts at e~^. The residue of the pole in 



Eq. (5.8) is 



It is this pole that is subtracted in the MS-scheme by adding to the Hamiltonian a counter-term 

b 



(5.9) 



(5.10) 



Similarly, the divergence arising from the contraction of two dipoles to a single dipole is subtracted by a 
counter-term 



,2 2e 
jJb 



II- 



with 



/ / 

J\x\<L J\i 



\v\<L 



Reducing this integral counter-term to a number, we subtract the residue of the single pole of 



\x\<L ly\<L 



■)-l I 



.kl +\y\ ) 



(5.11) 
(5.12) 

(5.13) 



la:l<L 1h|<Z, 



Note that the regulator L cuts off both integrations. One can now either utilize some simple algebra or show by 
the methods of conformal mapping (see [7]) that the residue is obtained by fixing one distance to 1 and by freely 
integrating over the remaining one 



Jo X J ~2-Z)p^_2D.j 



(5.14) 



As a result, the model is UV-finite at 1-loop order, if we use in the renormalized Hamiltonian (5.1) the renormal- 
ization factors Z and 



Z = l-{2-D) 



, I + Oib^) . 

e S 



Note that due to Eq. (4.19) no counter-term for 2 is necessary. 

The renormalized field and coupUng are re-expressed in terms of their bare counterparts through 

ro(a;) = Z^''^ r{x) , bo = bZb Z'^''^ . 



(5.15) 
(5.16) 

(5.17) 
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Finally, the renormalization group functions are obtained from the variation of the coupling constant and the field 

with respect to the renormalization scale ji, keeping the bare coupling fixed. The flow of the coupUng is written 
in terms of Z and with the diagrams given in Eqs. (5.9) and (5.14) as 



= -eh + 



-eb 



l + 6^1nZ6 + f6^1nZ 

- vd 



+ ) ) fe^ + 0(&3) 



(5.18) 



Similarly, the full dimension of the field (the exponent entering into the correlation function) is obtained as 



lh\ 2-g 1 d 



bo 



2-D 



1-6 



1 + 6 



2D 



(5.19) 



Note that minimal subtraction is used on the level of counter-terms or equivalently .^-factors. Since Z enters as 
Z"* into the /3-function, the latter also contains a factor of d in the 1-loop approximation, i.e. Z'' is not minimally 
renormalized. In order to calculate the leading order in e, the factor of d can be replaced by dc = ■ 

The /3-function has a non-trivial fixed-point with /3(6*) — 0, which has positive slope and thus describes the 
behavior of the model at large distances. The anomalous dimension i^* := i^{b*) becomes to first order in e 



2-D 



1 + 



2D 



1 ^{^r 



+ 1 



(5.20) 



For polymers, this result reduces to the well-known formula i^*(D = 1) = | + + 0((4 — d)^). 



6 Results for self-avoiding membranes from 2-loop calculations 

Two loop results are obtained by an explicit analytic and numeric integration of the following combination of 
diagrams (which are all of order 1/e) 




(6.1) 



One then has the freedom to extrapolate about any point on the critical curve e(D, d) = 0. This freedom - or 
ambiguity - can be used to optimize the results. The results of such a procedure for the 2-loop cacluations of 
V* are given on figure 6.1 for membranes (D = 2) in d dimensions (2 < d < 20). We see that for d —f oo 
the prediction of a Gaussian variational method Vva-r = ^ becomes exact. For small d, the prediction made 
by Flory's argument I'Fiory = 3^ is close to our results. This is a non-trivial result, since the membrane case 
corresponds to e = 4 and in comparison with polymers in d = 3, where s = 1/2, the 2-loop corrections were 
expected to be large. In fact, they are small when one expands around the critical curve e = for an adequate 
range of D ~ 1.5 (depending slightly on d and on the choice of variables) and a suitable choice of extrapolation 
variables. In this case the 2-loop corrections are even smaller than the 1-loop corrections and allow for more 
reUable extrapolations to g = 4. This can be understood from the large order behavior [10]. 
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Figure 6.1: Extrapolation of the 2-loop results in d and e for membranes D = 2md dimensions, using the expansion of {d + 2) (squares). 
The soUd line is the prediction made by Hory's approximation, the dashed line by the variational ansatz. 



Let us now turn to the physically relevant case of membranes in three dimensions (D — 2, d — 3). Our 
calculations predict an exponent « 0.85 or equivalently a fractal dimension of df « 2.4, which is in agreement 
with those experiments and simulations which find a fractal phase. However, this is still under debate, and a lot 
of evidence has been collected that the flat phase is generic. (For a more detailed discussion of this problem, see 
the review [7].) 

7 Outlook 

In these lecture notes, we have only been able to discuss the simplest applications of non-local field theories. 
Much more has been achieved during the last years. 

First of all, it has been demonstrated that the dynamics, described by a Langevin-equation 

^''^'''*) = ~^M)+^('^'*)' {a^,t)C{x',t'))=X6''ix~x')5{t-t') (7.1) 

also leads to a renormalizable field theory [11], where the dynamic exponent z, defined by the decay of the auto- 
correlation function ([r(a;, i) - r(a;, i')]^) ~ \t-t'\'^^^ isgivenby « = 2-|-Z)/i/* = 2-|-dt, where d/ = D/u* 
is the fractal dimension of the membrane. This result had been suggested long-time ago by DeGennes [12] for 
polymers and by Kardar et al. for membranes [13], but only with the methods discussed above, a proof of that 
conjectures could be given [11]. When hydrodynamics is included, the dynamical exponent z changes to 2 = d 
[14]. 

Interestingly, also disorder can be treated via the same methods, since averaging over disorder leads to in- 
teractions very similar to self-avoidance. The above methods have been applied to the motion of a polymer or 
membrane in a static disordered force field with both potential and non-conserved parts - leading to new universal 
physics [15,16]. 

Also the advection of a polymer in a turbulent flow has been analyzed [17], paralling the discusion of the 
passive advection of particles, which in the turbulence community is known as the passive scalar problem. 

Another interesting generalization is to anisotropically tethered membranes [18-20]. In these models, the 
membrane is more rigid in one direction, forming a tubular phase. 

One of the most useful tricks for self-avoiding polymers is the mapping onto a massive scalar field-theory, 
i.e. a 0* -theory [21]. In the limit, where the number of compoents n — » 0, resiflts for self-avoiding polymers are 
obtained. Stated differently, 0* -theory is a generaUzation of self-avoiding polymers. Another generalization has 
been discussed above: The generalization to membranes with internal dimension D ^ 1. The question arizes, 
whether a common generalization of both the oi^-theory and polymerized membranes is possible. Such a model 
has indeed been constructed [22,23], and leads to interesting new physics. 
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Another still open question is tiie analysis of the spectrum of subdominant operators [24], This may give a 
hint of why polymerized membranes are generically seen flat in simulations. 
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